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The expectation value of the complex phase factor of the fermion determinant is computed in the 
microscopic domain of QCD at nonzero chemical potential. We find that the average phase factor 
is non-vanishing below a critical value of the chemical potential equal to the half the pion mass 
and vanishes exponentially in the volume for larger values of the chemical potential. This holds for 
QCD with dynamical quarks as well as for quenched and phase quenched QCD. The average phase 
factor has an essential singularity for zero chemical potential and cannot be obtained by analytic 
continuation from imaginary chemical potential or by means of a Taylor expansion. The leading 
order correction in the p-expansion of the chiral Lagrangian is calculated as well. 
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I. INTRODUCTION 

A nonzero baryon density in QCD is achieved by introducing a chemical potential which enhances the propagation 
of quarks in the forward time direction and inhibits their propagation in the backward time direction. This imbalance 
makes the fermion determinant in the Euclidean formulation of QCD complex. The integration measure in the 
QCD partition function therefore cannot be used directly to define a probabilistic measure to generate gauge field 
configurations. Because of this sign problem, standard lattice QCD does not apply at nonzero baryon chemical 
potential. This does not mean that lattice QCD simulations at nonzero chemical potential are all altogether impossible. 
If the complex phase factor of the fermion determinant does not fluctuate strongly the sign problem may become 
manageable. The reweighting method [l| deals with the sign problem by including the phase factor in the observable 
rather than in the measure. While this method works in principle, its limitations are set by the strength of the 
fluctuations of the phase of the fermion determinant. It is therefore of considerable interest to understand how the 
fluctuations of the phase factor depend on the temperature, volume, chemical potential, quark mass as well as the 
lattice cutoff. 

In this paper we study the severity of the sign problem in QCD as expressed through the vacuum expectation value 
of the phase factor 

^2^e ^ det(D + hjq + m) 

det{D^ + ^i-fo + m)' ^ ' 

Since the average phase factor is the ratio of two partition functions, it is typically exponentially small in the volume 
requiring that the number of configurations needed in a reweighting approach is exponentially large with the volume. 

We will consider low temperatures where the chiral condensate E is nonzero and pions dominate the excitations of 
the QCD vacuum. We will take the square of the chemical potential inversely proportional to the volume 

M^F^y ^ 1, (2) 

where F^r is the pion decay constant. If furthermore 

toEF - 1 and FA^cD > 1, (3) 

the QCD partition function reduces to the zero momentum limit of the corresponding chiral Lagrangian. The scaling 
of fi and m in ([2]) and ([3]) defines the microscopic limit of QCD and is also referred to as the e-regime [1, Q . In this 
limit, the zero momentum modes of the pions factorize from the QCD(-like) partition function and exact analytical 
expressions can be obtained. We will derive exact expressions for the average phase factor in both quenched QCD and 
QCD with dynamical quarks as well as in its phase quenched version. The result in all cases is that the expectation 
value of the phase factor is exponentially small in the volume when /i > m^/2. For chemical potentials smaller 
than half the pion mass the average phase factor is nonzero. This suggests that the exact analytical predictions for 
/i < m7r/2 can be tested by means of lattice QCD even in the presence of dynamical quarks. 

In the past few years several lattice QCD methods have been developed to study the properties of QCD at nonzero 
chemical potential. In particular, progress has been made with reweighting 0, 0], irnaginary chemical potential 

I, B i [1] the Taylor expansion method P, [13, [HI [I2l ^^id the density of states method ^Ts^ (for a recent review see 
14|). Generally, one expects these methods to work if the average phase of the quark determinant fluctuates only 

mildly. For this reason several lattice measurements of the phase have appeared [H, [l^ [l^, [13, [3, [H, [13] ■ In this 
paper we will find that the quenched as well as the unquenched average of the phase factor has an essential singularity 
at zero chemical potential and cannot be obtained by analytic continuation from imaginary chemical potential or by 
means of a Taylor expansion. We give the explicit form of the non-analytic part and show that it is responsible for 
the abrupt change of the average phase factor at /x = m^/2. We also give the exact predictions for purely imaginary 

II. Since these are analytic in /i they can be tested directly with simulations at imaginary fj, or by means of Taylor 
expansion. We hope that the predictions for the average phase both for /i^ > and for jj? < Q will motivate lattice 
studies of the phase fiuctuations in the microscopic domain and thereby deepen our understanding of QCD at nonzero 
chemical potential. 

The scale u = m7r/2 has been part of lattice QCD at nonzero chemical potential starting from the very first 
calculations [13, [2l|. Although it manifests itself in different ways in quenched and unquenched QCD, it has the same 
physical origin, namely condensation of charged pions. 

The most direct effect of pion condensation occurs in QCD with an equal number of quarks and conjugate quarks 
(referred to as phase quenched QCD). For a conjugate quark, the fermion determinant is replaced by its complex 
conjugate, and therefore this theory corresponds to ignoring the phase factor of the fermion determinant. Since 
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complex conjugation of a fermion determinant corresponds to changing the sign of the chemical potential [24I, two- 
flavor phase quenched QCD is identical to QCD with a nonzero isospin chemical potential. Because the pions have 
a nonzero isospin charge and are weakly interacting they will Bose condense for fi > at zero temperature 

The quenched approximation at nonzero chemical potential is the zero flavor limit of QCD with an equal number 
of quarks and conjugate quarks [13, H^. It therefore inherits the scale /i = For example, one finds that the 

quark number density is nonzero for /j, > m7r/2 pTI. [29|. 

In QCD with dynamical quarks the Dirac spectrum for typical gauge field configurations has the same support as 
in the quenched case. Therefore, for quark masses such that m-^jl < fj,, one or more of the quark masses are inside 
the domain of the Dirac eigenvalues, and the sign problem is severe. It has to be like that because the phase factor 
has to wipe out the Bose condensate of the phase quenched theory (see also the review [M])- In QCD with dynamical 
quarks, the scale /i — is thus the scale where phase fluctuations become strong. This is also illustrated by the 

spectral density of full QCD: When the quark mass is inside the support of the eigenvalues the unquenched spectral 
density becomes a complex function with strong oscillations in a finite domain of its support [sil . [3^ . Physically these 
oscillations, with an amplitude that diverges exponentially with the volume and a period inversely proportional to 
the volume, are important because they are responsible for the discontinuity in the chiral condensate [33| . 

The difference between the phase quenched theory (where the phase of the fermion determinant is ignored) and 
unquenched QCD gives an indication of the importance of the phase factor. In fact the ratio of the unquenched 
partition function and the phase quenched partition function for two fiavors is exactly the expectation value of the 
phase factor in the phase quenched theory ([1]). For /i < m^/2 and zero temperature the free energy of the two partition 
functions is the same in the thermodynamic limit, but for fi > m^/2, because of the formation of a Bose condensate 
in the phase quenched theory, the two free energies become different. This results in an exponential suppression of 
the average phase factor. In this paper we will show that the same behavior is found for the expectation value of the 
phase factor with respect to the full QCD partition function and the quenched QCD partition function. We will do 
this by an explicit calculation of the average phase factor in the microscopic domain of QCD. 

In the microscopic limit QCD is equivalent to a chiral random matrix model with the same global symmetries 
[3^ [U, [3^, [3^ . The reason is that the static part of the chiral Lagrangian is uniquely determined by global symmetries. 
We exploit this equivalence to calculate the average phase factor in this domain. Alternatively, we could have started 
from the static part of a chiral Lagrangian and performed the integral over the zero momentum modes of the Goldstone 
fields. Since the phase factor includes a bosonic determinant at nonzero chemical potential this method is, however, 
technically demanding and has not been worked out up to now. 

In the microscopic limit it is useful to consider a fixed topological charge sector rather than a fixed vacuum angle 
Q . For simplicity, in this paper we will only consider the sector of zero topological charge. The results for arbitrary 
topological charge follow by an immediate generalization. 

Part of the results derived in this paper were announced in the letter |37|. In addition to providing details of 
the derivations, we obtain below explicit expressions for the average phase factor and show that it is a non-analytic 
function of the chemical potential. We also compute the corrections to the microscopic results to leading order in the 
p-expansion of the chiral Lagrangian. 

The behavior of the average phase factor can be understood from general arguments presented in the next section. 
The connection between these general arguments and the exact expressions in the microscopic limit are explained in 
section IIIII where we analyze the average phase factor in the phase quenched theory. In sections IIVI and |V] we derive 
exact expressions for the average phase factor in the microscopic limit. The thermodynamic limit of these results is 
evaluated in section IVH and the results obtained from general arguments will be confirmed. Section IVIIl contains a 
discussion of the analytic properties of the average phase factor, and we finish with concluding remarks. Technical 
details are worked out in two appendices. 

II. THE AVERAGE PHASE FACTOR FROM A GENERAL ARGUMENT 

In this section we will consider the phase fiuctuations of the fermion determinant in the QCD partition function 

Znj = {dct^f {D + Ai7o + ™) ) , (4) 

using mean field arguments at low energy. (Throughout this manuscript (. . .) denotes the quenched average. Dy- 
namical averages are labeled by a subscript referring to the number of flavors Nf.) In order to quantify the phase 
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fluctuations of the fermion determinant we compute the average phase factor |54l | 

\ aet(L'T + fj,^Q + m) / 

Notice that the average phase factor is the ratio of two partition functions: the partition function with an additional 
fermionic quark and an additional conjugate bosonic quark and the ordinary QCD partition function 

(e )Nf - — ^ • (b) 

This implies that the average phase factor is necessarily real. We will consider the low temperature limit with quark 
masses and chemical potentials well below Aqcd- Then the partition function can be described in terms of a chiral 
Lagrangian. In this section we will focus on the mean field limit where the fields are replaced by their saddle point 
values. In this limit the partition functions assumes the general form 



^ n 



1 



k 



(7) 



where the Jacobian, J, is from the measure of the Goldstonc manifold at the saddle point, and mT^(fi) are the /i 
dependent masses of the Goldstonc modes. (For a discussion of the parameterization dependence of the masses of the 
Goldstone modes see Appendix [XI) The free energy density at the saddle point, F, is intensive, and the prefactor can 
be written as a 1/V correction to the free energy. 

The QCD partition function, Z^f, in this approximation is completely independent of the chemical potential; since 
the pions have zero baryon charge neither the free energy density nor the exponential prefactor depend on fi. 

At zero chemical potential the fermion determinant is real so in this case the two additional determinants in 
ZNf+i\i* cancel and the average phase factor is 1. When the chemical potential is nonzero, the two additional flavors 
become important. The combinations of the conjugate bosonic quark with any of the fermionic quarks forms a charged 
Goldstino. Because of the condensation of charged Goldstinos for /i > m7r/2 we need to consider separately the cases 
/i < m7r/2 and /i > 711^^/2. 

For /I < m7r/2 the free energy F is independent of the chemical potential. In the expectation value of the phase 
factor the contributions of the free energy to the average phase cancel. The prefactor in ([7]) however depends on /i 
because the /^-dependence of the masses of the charged Goldstone modes. Half of the 2(iV/ + 1) charged modes have 
masses TO,r(/i) = m,r + 2/^ while the other half have masses rn^(/i) — TO^r ~ 2/^. The mean field result for the ratio of 
Zni+i\1' and Znj is thus given by 

(e = 2iv7+2 = (1 - ;;;^) ' fo^' ^ < o'"- (8) 

For 2/i > TOjr the free energy of 2^Af/+i|i* depends of /i. The free energy has the usual /i dependence for the static 
part of a chiral Lagrangian in the condensed phase (see for example (23l [24. l26t ) . Subtracting the free energy at /i = 
we find 

(e^^^)^, ~ e-2^^^^^^(i-™^/4'^^)^ for ^, > \m^. (9) 

The average phase factor is thus suppressed exponentially with the volume for ^ > m-n/l. Care has to be taken to 
compute the exponential prefactor. For > nii^jl there are massless modes associated with the condensation of pions. 
Consequently the leading contribution to the prefactor given in ([7]) vanishes. One way to compute the subleading 
contributions to the prefactor is from the thermodynamic limit of the exact microscopic result. In the next section 
we perform this calculation in the phase quenched theory. 

In section IVII we will confirm the results obtained here from the thermodynamic limit of the exact microscopic 
expressions. 



III. THE AVERAGE PHASE FACTOR IN THE PHASE QUENCHED THEORY 



In this section we will discuss the chemical potential dependence of the average phase factor in the phase quenched 
theory, which can be easily obtained from known results. First we will discuss the exact microscopic result and then 
relate it to the mean field results of section [III 
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The phase quenched partition function with two dynamical flavors is just the two-flavor QCD partition function 
without the phase of the fermion determinant, 



^1+1 



-(m;/i) = (I det(i5 + ^70 + m)|2). (10) 



This partition function is equivalent to QCD at nonzero isospin chemical [24]. The expectation value of the phase 
factor of the fermion determinant for the phase quenched theory is thus given by 

I 2^e^ (dct^(i:' + ^7o+m)) ^ Zn;=2 . 

^ (|dct(I? + Ai7o + m)|2) ■ ^ ' 

Its microscopic limit therefore follows immediately from the microscopic limit of the two partition functions. 
The microscopic limit of the phase quenched partition function is given by [38| 

/■I 

Zi+i.(m;/i) = 26^^' / dtte-'^l'^^'' ll{mt), (12) 
Jo 

where 

rh — mVY, and fl — fiF^rVv. (13) 

The two-flavor partition function is /i-independent at scales well below the nucleon mass. Its microscopic limit is 
given by 3 

ZNfMm) = llirh) - Ii{m). (14) 
The /i-dependence of the phase factor thus follows from the /i-dependence of the phase quenched partition function 

m 



We now take the thermodynamic limit, m — > cx) and /i^ ^ 00, of the microscopic result for the phase factor in the 
phase quenched theory. In this limit the phase quenched theory has a second order phase transition at rh/{2ji^) = 1. 
To show this we calculate the integral in (|12p by a saddle point approximation. The transition between the two phases 
takes place when the saddle point hits the boundary of the integration region, i.e. when 

If we use the Gell-Mann— Oakcs— Renner relation we find that this corresponds to the expected critical value of the 
chemical potential oi — m^/2. 

For < m^/2, the saddle point is outside the integration region and the leading contribution to the integral comes 
from the edge of the integral at t = 1. In the thermodynamic limit we obtain (this is derived under the assumption 
that the integral is cut off by the exponential factor rather than the Gaussian factor, which in the thermodynamic 
limit is violated only very close to the critical point) 

Z.,H-;A)-^-(;^^e- for .n>2A^ (17) 

We find that for m > 2//^ the free energy of the phase quenched theory is independent in the thermodynamic limit. 
The thermodynamic limit of the microscopic two flavor partition function is 

Z^,=2(m;A)-^e2'^. (18) 

In the average phase (jlSp the free energies cancel. The /x-dependence only enters through 1/V corrections, i.e. through 
the pre-exponential factors 

(e^^«)i+i.^(l-¥)e° = (l-l^)eO, (19) 
m mi. 
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where we made use of the Gell-Mann— Oakes— Renner relation. This is precisely what we find from the general 
argument given in previous section. To see this we evaluate the two partition functions using ([7|). First, note that the 
particle content of the phase quenched partition function is a charged Goldstonc boson and its anti-particle and two 
neutral Goldstone bosons. For /j, < m.n./2, their masses are given by m^r — m^r + 2/i, m^r, m-^, respectively, and 
the free energy is equal to by 2mE. The Jacobian in the integration measure is a constant. The ordinary two flavor 
partition function is of course obtained by setting ^ — in the expression obtained for the phase quenched partition 
function. 

For > m,r/2 the thermodynamic limit of ^2]) allows us to determine both the free energy and the exponential 
prefactor. In this case the saddle point is inside the integration region and the saddle point approximation to the 
partition function (jl2p is given by 

Zi+i.(m;A) ~ \/^-:Ve"'/2A^+2'^' for rh < 2fx\ (20) 
V ZTT 



Notice that the prefactor of the exponential is larger by a factor of VV than for fj, < mTr/2. The reason is that the 
Goldstone mode with mass — 2/i becomes exactly massless for 2/i > m^. To use the saddle point approximation 
one must introduce a pion condensate source term which lifts this mass away from zero (see Appendix |X| . 

The two flavor partition function remains fi independent so that in the thermodynamic limit the average phase 
factor is given by 

f (1 - 1^ ) e° for > 2fi, 

(e^'^)i+i. ^ <^ " (21) 

1 2MF^e-2VF=p^(l-mJ/V)^ for TO^ < 2/i. 

Below we will derive the exact microscopic expression for the average phase factor in the quenched and unquenched 
theories. To perform the calculations we will make use of recently developed random matrix techniques. 



IV. THE RANDOM MATRIX MODEL 



In order to compute the average phase factor in the microscopic limit we use a random matrix approach. In this 
limit the mass and chemical potential dependence of the QCD partition function is given by a random matrix partition 
function (see [s^ for a review). The expectation value of the phase factor is then given by the ratio of two partition 
functions 

/ 2»0\ ^jV/+i|i-(™;A^) . 

(e )n, = — 7} — -, ^ — : (22 

where, in general, the random matrix partition function with Nf quark flavors of mass my and one pair of a regular 
quark and a conjugate bosonic quark both with mass m is defined by 

Nf 

^iv,+i|i.({™/},"^;A^) ^ |d<&d*^^($)^„^(vI/)JJdet(P(/.)+my) (23) 

In the computation of the average phase factor we always take the quark masses to be degenerate. The Dirac operator 
is given by 

^<'')^(.*.+°,* ^ (^^> 

Here, $ and ^' are complex N x N matrices with the same Gaussian weight function 

wciX) = exp(-A^TrXtX) . (25) 

The partition function Zn^ (m; fj.) with Nf flavors of mass m is given by (j23p without the ratio of the two determi- 
nants multiplying the Nf flavor quark determinant. The microscopic limit of the random matrix model is obtained 
by taking the size of the matrices infinitely large, N oo, while keeping the products 

m = 2mN and /i^ = 2fi^N (26) 
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fixed. Expressions involving the microscopic partition function are denoted by explicitly writing the microscopic 
variables fi and m. 

The random matrix partition function (|23p can be rewritten in terms of an integral over the eigenvalues of T) and 
the unitary matrices that diagonalize T). The integral over the unitary matrices can be performed analytically [3^ 
resulting in the eigenvalue representation of the partition function 



N 

ZNi+i\i*{.{mf},m;^j.) ^ j \\_d?ZkVNi+i\i*{{z,},{z*};ii), 
■^^ k=l 



(27) 



where the integration extends over the full complex plane. The joint probability distribution of the eigenvalues is 
given by 



V 



Nf + l\l' 



^jv \^Ni{zf})\ Y[ w{zk,zl;ii) W_{m) - zl) 
^ fc=i /=i 



("^^ - zl) 



The Vandermonde determinant is defined as 



N 



/^M{{zf}) ^ n - 



and the weight function reads [32 



The eigenvalue representation makes it possible to define orthogonal polynomials in the complex plane 
I42I . l43i . i44i . ,45j . In order to evaluate the average phase factor we will make use of such orthogonal polynomials and 
their Cauchy transform. The complex Laguerre polynomials given by [32| 



(Zk + Zk* ) 



(28) 



(29) 



(30) 



Pk{z;n) 



N 



k\Lk 



Nz^ 



(31) 



are the orthogonal polynomials corresponding to the weight w{z^ z*; /i) given in (j30p . They satisfy the orthogonality 
relation 



(fzw{z,z*;^)pk{z;fj,)pi{z;fi)* = rk 6ki 



with the norm 



rk 



7r/^2 (l+^2)2fc ^, 



The Cauchy transform of the orthogonal polynomials is defined as 

1 



d" z- 



w{z,z*;^)pI{z;^), 



(32) 



(33) 



(34) 



where C indicates that the integration extends over the complex plane. 



THE AVERAGE PHASE FACTOR IN THE MICROSCOPIC LIMIT 



Partition functions that are given by expectation values of ratios of determinants can be expressed simply in terms 
of orthogonal polynomials and their Cauchy transforms [Sll, |4l|, \^ |4^ 13, ■ Generalizing the results of [i^ we 
find that the partition function ()23p for matrix size N is given by 



-■Nf + l\l 



.({m/},m;/i) 



rN^iAN^.+i{{mj}) 



hN-i{m;fi) 

PN-i{m) 

PN-i{mi) 



PN+Nf {m) 
PN+Nf{mi) 



PN~i{mNf) ■■■ PN+NfirriNf) 



(35) 
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Using the orthogonality of the polynomials on which this integral is based, it can be shown that this integral vanishes 
El 



xiJ)(m;A) = 0. 



(44) 



What remains is the 5-function part of (H^ in combination with the part of the of factors in the integrand that is 
continuous accross the negative imaginary y-axis. By inspection one easily finds that the contribution of the pole at 
s = —X — m vanishes. For the contribution of the pole at s = a; — m we have to distinguish rh — x < x and m — x > x. 



Because v^*/ has a cut for all values of s (notice that z* — a; — s =p e), in the first case the nonvanishing 
contribution comes from the cut in Xodzp) and in the second case from the cut in Kq^—z*). Wc thus find 



xW(m;A) 



4/i2 



re e 



™/2 ((2x-m)m 
ax(2x — m)lo 



m/2 



dx(2x — rh)Ki 



( {2x — rh)m 
[ 4A^ 



exp 



4^2 
(2x - to)2 



exp 



im-2x)^\ 



Si_,AKo{m-2x)] 



(45) 



This result can be simplified to 



4A 



duueyi^[-^]Ko ( ^ ) (udufloiu) + / duu exp[--^]/o ( ^ ) {udu)''Ko{u) 



.(46) 



Making use of the identities derived in Appendix [B] these integrals can be rewritten as the sum of a polynomial in /i^ 
and an integral with an essential singularity at /i = 0, 



1 



X(")(m;/i) = Ko{m) - e-^^^ 

X(i)(m; /i) = -mKi{m) + ifi^Koim) - e"^''"-^e 



duuexp[--—^]Ko ( — I Io{u) 



1 



4/t2 



duuexp[-—]Ko — — I uduIo{u) 



8/i 



4/i2 



-e^^^^ -—^e~^ / duuexp[-^]/s:o ( ^ ) {udufh[u) 



4^2 



8A2 



4^2 



X(3)(to; /t) = Ka{rh){l2rn^ fi^ + 64/i^ + + 2m2 + l6/i2) _ Ki{rh){A8rnji^ + + 24m/i2) 



1 _j 



4^2 



8^2 



4/i2 



(47) 



(48) 



(49) 



(50) 



The analytic properties of the X'^^^ will be discussed further in section IVlTl 

In oder to complete the computation of the average phase factor (|22p we only need to recall that the ordinary (/i 
independent) microscopic partition function is given by ^3] 

Z^,(m) ~ m-^/(^/-i) det[4+'/o(m)]fe,=o,...,7V,-i. (51) 
In the next section we will discuss explicitly the result for the average phase factor for Nf = 0, 1 and 2. 



VI. EXPLICIT RESULTS AND THE THERMODYNAMIC LIMIT OF THE PHASE FACTOR 

In this section we will take a closer look at several specific cases and derive the large m and large jl asymptotic 
expansions of the average phase factor from the exact microscopic expressions given in the previous section. We 
refer to this as the thermodynamic limit of the microscopic results. These results confirm the expressions for the 
average phase factor obtained in section |TT] from the mean field argument. For the technical details of the asymptotic 
expansion we refer to Appendix [Ul 
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1. The Quenched Case 
We first consider the expectation value of the phase factor in the quenched case. Then 



[e )Nf=0 



In{in) mli{rh) 
l-ifi^Io{m)Ko{m) 



—-e ^ / da;xexp[-— ]Xo — ^ I (/o(x)to/i(to) - a;/i(x)/o(m)) , 



(52) 



(53) 



where we made use of the Wronskian identity rhloijh) Ki{m) + jhli{rh) Koirh) = 1. 

We now take the thermodynamic hmit m —^ oo and fi —^ oo. In the normal phase we can use the leading order 
asymptotic expansions for the X^^^ and the Bessel functions whereas in the condensed phase we have to include the 
subleading corrections. Substituting the expressions obtained in Appendix [Cl we find (the result for m-^ < 2fj, is not 
valid very close to the critical point (see Appendix C)) 



(1 - ^) e 



[e )Nf=o 



1 1 1 



' for m-^ > 2fi, 

-2yF,=M=(i-m^/V)= for < 2^l. 



(54) 



For rn^ < 2fi the leading order terms contributing to the prefactor cancel. This asymptotic result agrees with the 
result we have obtained from the general arguments in section HIl 



2. Full QCD for Nf = 1 

For one flavor, the explicit microscopic expression for the expectation value of the phase factor is given by 

1 



(e )Nf=i 



2to^/o(t 



loifh) mli{m) m^/o(m) 

mli(rh) m^/o(m) ?7i^(2/o(m) + rh/i(m)) 



(55) 



where the AT'*^' have been given in the previous section and we used that Zj\[^^i{fh; fi) — Io{rh). With the expressions 
found for A^*^) the average phase factor can be written as 



him) 



{e'nNf=i = 1-4 mKoim)Io{m) + -^l-l-{Ko{m) - mK^{rn)) f/ + 8Ko{m) \ , . ' ft 



dxx exp[-^^]A:o 



Io{x) 5.Ji{x) (5^fh{x) 
Io{m) 5raIo{m) {6rhYlo{rn) 



(56) 



The large m and // limit of the expectation value of the phase follows from the asymptotic expressions for the A^'^^ 
given in Appendix [C] and the asymptotic expressions for the Bessel functions 



(1 



for > 2/i 



(e )jV/=i 



(57) 



0{-^)e-^VF-^.^il~rnl/4^.-f for < 2^ 



The first three orders contributing to the prefactor of the result for < 2/i vanish. 

Since the average phase factor is the ratio of two partition functions it is necessarily real even if the "statistical 
weight" is complex such as for Nf — 1. This does not imply that (sin2^?) vanishes such as in the quenched theory or 
the phase quenched theory. Since (sin 20) — {exp{2i9) — exp{~2i6)) /2i and 



, 2^e^ ZNf = l'{m-^l) 

{e )Nf=i = — 7 r = -•-) 



(58) 
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we have in the microscopic hmit that 

(sin 20)Ar,=i = ^ [{e''')Nf=i - 1] . (59) 

We thus find that the expectation value {sm29) f^^^i is generahy nonzero and purely imaginary. Its thermodynamic 
limit is given by 

f ^((l-^)'-l) for m^>2^,, 
{sin{2e))N,=i ^ I " (60) 

[ for < 2fi. 

Since the average of the inverse phase factor is 1 we automatically also know the variance of the phase factor 

\e e )ns=\ — \& )Nf=i\e )Nf=i — ^—\e )Nf=i- (olj 

So while the average phase factor goes to zero, its variance goes to 1 for /i > 771-^/2, what would also have been the 
result for a uniformly random distribution of the phase. This central feature of the sign problem is also present for 
larger number of flavors as well as for quenched and phase quenched QCD (in the latter two cases there is of course 
no sign problem). 



3. Two dynamical flavors 



With two flavors the explicit microscopic expression for the average phase factor is 

\e )Nf=2 



8m6(/o(m)2 - /i(to)2) 

X(o)(m*;A) XW(m*;fi) X(^\m*-p) X(3)(m*;A) 

Io{m) Sfnloim) (5rhfIo{m) {Srh)^Io{'ih) 

Srhloim) {5rhflo{m) (5rhfIo{m) {6rh)^Io{in) 

{5rafh{rh) (Sfnfloim) (<5™)%(m) (<5™)5/o(m) 



(62) 



For fi < m-^/2 the thermodynamic limit of this result is 



(1 



for > 2/z. 



(63) 



As the chemical potential increases beyond half the pion mass the average phase factor is again exponentially small 
in the volume in agreement with The fast convergence to the thermodynamic limit is illustrated in figure [TJ 



VII. ANALYTIC CONTINUATION FROM IMAGINARY CHEMICAL POTENTIAL 

For purely imaginary chemical potential the fermion determinant is real, and the phase factor is equal to unity. 
Thus it seems to be impossible to calculate the average phase factor at imaginary values of fi and then perform an 
analytic continuation in fi. However, if we first use that det(Z? + /_t7o + m)* — det(Z? — fi"fQ + m) to write the average 
phase factor at real chemical potential as 

g2»e(p)\ _ / ..^H^ -r ^7U -r v ^ ^^^^ 



\ / Nf \det{D - fi'jQ + m) / 



and subsequently substitute fi by ifi we obtain a real number 



\ det(i:> - i^7o + m) / 



N, 



that is (typically) different from one. Below we will use this analytical continuation between real and imaginary 
chemical potential. 
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FIG. 1; Left: The average phase factor for two dynamical flavors as a function of the chemical potential for fixed quark mass. 
The full curve represents the result in the thermodynamic limit. Note that the convergence to the thermodynamic limit is 
particular rapid for small values of the chemical potential. Right: The average phase factor as a function of the quark mass 
for fixed chemical potential. The full curve again displays the result in the thermodynamic limit. 



The microscopic limit of the quenched av erag e of ((65|) for imaginary chemical potential can be derived by means of 
the supersymmetric method. It is given by [5a | 



e 



)iv,=o = 1 + 4A2jo(^)i^o(m). (66) 



Since the result is a polynomial in /i^, analytical continuation to real chemical potential is trivial: we simply flip the 
sign of fj?. This result should be contrasted with the result for real /i given in Eq. ((53)) : 



i-co / xrh\ 

(e2*'(^))jv,=o = l-4fi^Io{m)Ko{m)-^^e ^ J dxxe ^ i^^^ \ {Ioix)mh{m) ~ xh{x)Io{m)) 

^ l-A^^Io{m)Ko{m) + il^Ko(^^^e'^'^^\loimf -h{mf) for f, ^ 0. (67) 

What is reproduced by the expression (j66p for purely imaginary chemical potential are the analytic terms in fj,, but 
the term with the essential singularity is not obtained. Note that the non-analytic term does not contribute to the 
Taylor series of (j67p . The Taylor expansion around /i = of (|67p is thus given by the first two terms only. This 
expansion can be analytically continued in and necessarily reproduces the result for imaginary /i because its Taylor 
expansion is finite. 

The microscopic limit of the average phase factor for imaginary fj, can also be evaluated for Nf ^ using the 
supersymmetric method of (48.] . However, it is possible to obtain {exp(2i9{i^)))Nf fi'om the results at real /i without 
a detailed computation. We only need that the average phase factor for imaginary fi is given by a finite power 
series in /i^. This immediately follows from the super symmetric approach where ifi only occurs in combination with 
Grassmann variables. Therefore, also for Nf ^ 0, the exact average phase factor at imaginary /i can be obtained from 
the analytical continuation of the Taylor expansion of the expression for real fi. For example, for Nf = 1 we obtain 
the result 

(e2'^(^^))Ar,.i = 1 + 4 (mKo{m)Io{m) + ^-0l{Ko{m) - mK,{m))) fi^ + ^K^{rh)^-^^^%^^^^ ^\ (68) 

Conversely, if we would have calculated the "phase factor" for imaginary chemical potential and then made the 
analytical continuation to real /i we only would have obtained the polynomial terms in fi. Terms for which all 
derivatives at ^ = vanish cannot possibly be obtained from a Taylor expansion. 

For ^ <^ tItt/S the term with the essential singularity is not important. However, it becomes the dominant term 
when /i approaches For /i > m-^jl the non-analytic term cancels the terms that are polynomial in /z, and is 

responsible for the exponential suppression of average phase factor. In the thermodynamic limit, the non-analytic 
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term has a phase transition at ^ = For ^ < 171^^/2 the saddle point is outside the integration range and the 

leading contribution to the integral comes from its lower limit. On easily shows that the integral vanishes for V oo 
and /i < For fj, > m,r/2 the saddle point is inside the integration domain and its leading order exactly cancels 

the thermodynamic limit of the terms that are polynomial in 

A similar essential singularity at = is found for Nf = 2, and the results obtained by Taylor expansion or analytic 
continuation in only reproduce the finite polynomial series in /i^. What is not reproduced by analytic continuation 
are precisely the contributions given by the integrals in (|T7)l - (|49p . 

The integrals with the essential singularity are due to Dirac eigenvalues with real part larger than m. This is best 
illustrated by considering X^"^ which is the partition function for one bosonic flavor [43|- In terms of a joint eigenvalue 
distribution it is given by 



If 2 ^ 1 

^ k fc=i 



2Y> (69) 



4y 



where the Vandermonde determinant and the joint probability distribution are defined in (|29p and pop . respectively. 
If we single out one eigenvalue, which we will denote by w this can be re-expressed as 

where we have used the identity (45j 



Af , N , N 



J. J. 7772 ^2 ^ ^ ™2 ^2 J, J, ^2 ^2 ('^ ) 



k=l 1^ j=l 3 k^j 3 



Using that in the microscopic limit the partition function Zj^^^^{w*) is given by 

Z^;l,{w*)~.I,{w*). (72) 
we obtain for the microscopic limit of Z^^^_-^ (which as before is denoted by X^^^ using the same normalization) 

The expression for X^^') given in (|47p is obtained by integration over the imaginary part of li by means of Cauchy's 
theorem and introducing the variable u = 2Re(w) — to. The integral over it > to in (j47p . thus corresponds to 
contributions with the real part of one of the eigenvalues larger than to. The same argument can be made for 
X('=>o)(to;A). 

For /t^ > to/2 the probability that one of the eigenvalues has a real part larger than to remains nonzero in the 
thermodynamic limit so that the non- analytic term becomes important. In the thermodynamic limit it exactly cancels 
the analytic term resulting in a vanishing average phase factor. In figure [2] we illustrate that the contribution of the 
non-analytic term (given by the difference between the upper and lower dashed curve) is important close to the critical 
value of fi and beyond. For comparison we also give the thermodynamic limit of the average phase factor (full curve). 

The source of the non-analyticity is the inverse determinant of a non-hermitian operator. For example, for the 
partition function with one bosonic flavor with imaginary chemical potential the commutator term in the chiral 
Lagrangian vanishes resulting in a microscopic partition function that is independent of /z. For real fi the bosonic 
determinant has to be regularized resulting in a /z-dependent effective Lagrangian (see (47| for details), and a partition 
function that is /z-dependent and therefore non-analytic at /i = 0. 

The message we would like to convey is that the Taylor expansion in ^ about zero of averages which involve an 
inverse quark determinant has zero radius of convergence. Also traces of the inverse of non hermitian operators lead 
to essential singularities at /x = unless the inverse is compensated by a factor from the fermionic determinant. The 
simplest example is the quenched chiral condensate which is non-analytic in fj, at fj, = 0. 

Note that the average phase factor in the phase quenched theory (llip does not involve averages of inverse determi- 
nants and is analytic in /i at = as can easily be verified from (|15p . 
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FIG. 2: The quenched average of the phase factor. Shown is the microscopic resuh for real /i (upper short dashed curve) and 
the microscopic result obtained by analytic continuation from imaginary chemical potential (lower long dashed curve) . The full 
curve represents the thermodynamic result for real fi. Note that while the average phase factor is non-analytic at p = 0, the 
analytic continuation from imaginary /i follows closely the correct result for 2/i/m,r "C 1. 



VIII. 1-LOOP CORRECTIONS 



So far we have focused on the average phase factor in the microscopic domain where m^F^ ^ l/V and fi^F^ ~ 1/V 
as the volume is taken to infinity. A perturbative expansion scheme where p ^ 1/i, TOttJ^tt ~ 1/-^^: ^ is 
known as the e-expansion The microscopic domain is the leading order term in the e-expansion. In this section 
we consider the so called p-expansion where p ^ 1/i, ^Jt-tt ^ 1/L, fi ^ l/L and work out the average phase factor in 
an expansion in m^L^ and ^^L^ for a finite but large box. In the p-expansion, the microscopic variables y?F^V and 
mTiV are large, and consequently the leading term in the expansion is given by the asymptotic limit of the microscopic 
expression. Below we will compute subleading terms up to order mj/i^i^ and ji^L^ . This corresponds to the 1-loop 
corrections generated by the nonzero momentum terms of the Goldstone bosons. For simplicity we will consider only 
the average phase factor in the phase quenched theory for /i < m7r/2. 



A. 1-loop Integral at ^ 7^ 

In the standard notation of chiral perturbation theory we have (ioj 

(e'^'>i+i- = = exp(Go(Ai = 0) - G^{^i)). (74) 

Zi + i* 

Notice that the contribution of the neutral Goldstone bosons cancels, and that we have a factor 2 in the exponent 
because the contribution of the two oppositely charged Goldstone bosons is the same. The 1-loop contributions of a 
single charged Goldstone boson with charge 2 (/i is the quark chemical potential so a pion made from an up and a 
down quark has charge 2) in a box with volume L"^ is given by 

e«°(^)/2 ^ exp[-i ^logbt + ml + {p,o 2*^)^)], (75) 

Pka 



where 



Pka = ka integer. (76) 
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(For a fermionic Goldstone boson the sign of the exponent is positive.) The contribution to the free energy given by 
G'o(/i) is divergent. However, the difference between Go{fj,) and Go{fi) for F = oo is finite, that is 

Go(/i) = Go(/i)|y=oo + .go(M) (77) 

with 5o(/^) finite. Since Go(/i)|v=oo is independent of /i for < (in dimensional regularization, see e.g.[50j) we 

have that 

(e2»«)i+i. ^ |X£z^ = exp(5o(M = 0) - goU^))- (78) 

Zl + i* 

We wish to evahiate this 1-loop contribution including the 1/V corrections. This can be done along the lines of [4§| 
where the expression for go was worked out for /i = 0. The first step is to perform a Poisson rcsummation 

Coin) = -J2^og{pl + ml + {pko~ 2z/i)2) ^-VJ2 f ^e^^"""'" log{f + ml + {po - 2iy.fl (79) 

where the sum is over all integers. The thermodynamic limit, Go|y=ooj is given by the term 1^ = 0. Since this term 
does not contribute to (exp(2i0)) we exclude it from now on. For /i < rriT, /2 we can shift the po variable by 2i/i. This 
results in 



e— 



e .In \^~'^ 



(80) 



where 7 is the Euler constant. The 1/e — 7-terms do not contribute because the term with all = has been excluded 
from the sum. After performing the momentum integrals, the limit e — > can be taken safely resulting in 

Next we change to dimensionless integration variables by 

A->A— with L = (nL^y/'^ (82) 
47r 

a 

and interchange the sums and the integral to arrive at 

roc J \ ,2 

go ^ ^A-''/2e-™Ji^A/4.(-Q^g-2M/oio5„Og--^ (83) 

^ a 1^ 

The integral over A can be split into into a part with A < 1 and a part with A > 1. On the second part we apply 
Jacobi's imaginary transformation 

oo I — oo 



^ e-''~'"^ = J-e^ J2 e-^-^. (84) 



a 

k— — oo 



This leads to 



+ r ^e^e-'"^^^^V4.(-Q^,-2../o^..o,-<^ _ 1) + (85) 

The Jacobi imaginary transformation is applied to the full sum so that the subtracted term before and after the 
transformation is different. This difference, denoted by tq, is given by 



ro 



A ./i A 
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Finally, we change A ^ 1/ A in the Jacobi transformed terms in go to obtain 



50 = 



e AI,^ — 1) 



^ dA 
T 



a I 



2 



(87) 



The expansion of tq follows immediately from [49| . The only modification is that the Goldstone boson mass in one 
of the terms now depends on the chemical potential. We find 



, {ml - 
ro = -fog 7 



4tt 



1 fmlL^'' ' ' 



^ — ^ 7i}n. 



{ml - 
47r 

2 t2 °o 



2 V 47r 

For the difference — tq we obtain the expansion 



^"2 



1 m^L 
- H ^ 

2 47r 



y L 



n\(n - 2) 

n=3 ^ ' 



mlL^ 
An 



90 - ro 



V V L 

Til ("On 



^ n!(2TO)! 

ra=0 m=0 ^ ' 



47r y ^ ^ ' 4^ 
For a d-dimensional hypercubic box the shape coefficients Sk,i are given by 



(89) 



Jo ^ J J 
S,,o - /'^A*^[(5:e--'?/^)'^-l]. 



(90) 



The shape coefficients Sk.o already enter in the expansion for /i = [49|. In the table below we give the numerical 
value of a few low order coefficients. The lowest order terms of the expansion are given by 



Sk,i 


/ = 


I = 1 


I = 2 


k^O 


0.09375685 


0.025131706 


0.025139306 


fc = -1 


0.11745759 


0.031290174 


0.031298318 


fc = -2 


0.15365038 


0.040637153 


0.040645920 


fc = -3 


0.21251023 


0.055744392 


0.055753879 


fc = -4 


0.31550608 


0.082023256 


0.082033581 



50 



(ml - 4Ai2)i2 



{ml - 4fi^)L^ 



4tt 

(ml - 4^2)^2 



7 - - + 5o.o + 5-2,0 + 2 
(l-5_i,o) + ^4^(l-5_i,o) 



1 ^mlL' 



m|L2 
87r 



47r 



47r 

(-4 + 2^-2,o)+( ^ 



47r 



:<5'o,o 



, m^L^ 3 



3,l) + (4^2^2)2(^^_^_^^^5_ 



4,2 - ^5_3,l) 



(91) 



The fi independent terms do not contribute to the average phase so that at next-to-next-to-leading order we find 

(92) 



e 1+1* 



4^2 /4^,^L^ ^ %ml^l^L^-lQ^,^L\ 11 



4 2 



.^4^^L2(-5.i.i + 5_3.l) - (4^2L^)2(^5_2,2 + ^5-4,2 - ^5-3.l) 
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FIG. 3: The average phase factor in the phase quenched theory as a function of the chemical potential in the normal phase. 
The full line gives the result at leading order (|2ip and the dashed line includes the corrections up to next-to-next-to-leading 
order (f92l) when m-KL — 1. 

The 0, wLtjL —>■ limit of this expression should reproduce the fi —> oo, m oo limit of the microscopic result 

which is indeed the case (see (|2T|) ). The correction factor computed in (|92p drives the average phase factor closer to 
1 (see figure [3]). Note that the correction is small even though m^L is set to 1. 

IX. CONCLUSIONS 

We have computed the average phase factor in the microscopic domain of quenched, unquenched as well as phase 
quenched QCD and have found that in the thermodynamic limit the average phase factor undergoes a phase transition 
at ji = For fj, < r7i^/2 it falls off polynomially with fj, and is equal to zero starting from fi — r7i^/2. This 

implies that the phase problem becomes severe for /x > 171^^/2. Physically this should be the case because the phase 
factor has to wipe out the Bose condensate that would be the vacuum state if the fermion determinant would have 
been replaced by its absolute value. 

The results for the average phase factor in the microscopic limit of QCD were derived from a chiral random matrix 
model exploiting the equivalence of these two theories in the microscopic domain. Powerful random matrix methods 
such as the method of orthogonal polynomials and their Cauchy transforms enabled us to obtain exact analytical 
results. The thermodynamic limit of these exact results, was also derived from simple mean field arguments. The 
starting point was that the average phase factor is the inverse ratio of the QCD partition function and the partition 
function with the fermion determinant of QCD and an additional conjugate bosonic quark and a fermionic quark. For 
fi < m^/2 the partition function with the additional quarks is in the same phase as the partition function without 
additional quarks. The ratio of the two partition functions is therefore determined by the 1/V corrections to the free 
energies. For larger values of /Lt the presence of the conjugate bosonic quark induces a pion condensate and the two 
free energies no longer cancel in the thermodynamic limit. This causes the exponential suppression with the volume 
of the average phase factor. 

While the average phase factor is exponentially small for /i > m^/2 the fluctuations of the phase factor becomes 
unity for fi > m^/2. The sign problem is therefore particularly severe for /i > to^/2. 

The finiteness of the average phase factor for ^ < suggests that it should be possible to test the predictions 

of this paper by means of lattice QCD simulations. We note however that the quenched and unquenched phase factor 
is non-analytic in /i at = 0. As has been shown explicitly in the quenched case, the phase factor as defined by the 
ratio of two determinants with opposite sign of the chemical potential, cannot be analytically continued from negative 

to positive jj,^. It is our experience that analytic! ty at /i = is lost for averages involving the inverse (determinant) 
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of the nonhermitian Dirac operator when the singularity is not compensated by a similar factor from the fermion 
determinant. In such cases lattice methods that rely on analyticity cannot be used. On the other hand, we expect 
observables that are derivatives of the usual fermionic partition functions to be are analytic in fj, for — > and can 
be computed by means of the Taylor expansion method or the imaginary chemical potential method. 

Since each eigenvalue contributes to the phase factor it is perhaps surprising that the average phase factor is not 
sensitive to the ultra-violet cutoff. The reason goes back to the renormalizablity of chiral perturbation theory: If the 
theory is renormalized at /i = 0, the fj, dependence does not introduce any additional infinities. Since the average 
phase factor is a partition function at /i divided by the same partition function /i = 0, it is therefore necessarily 
ultraviolet finite. We have shown this explicitly to 1-loop order for the phase quenched partition function. This 
follows from the microscopic theory as well. In lattice QCD the Hasenfratz-Karsch prescription to introduce the 
chemical potential does not give rise to additional ultraviolet divergences in lattice perturbation theory. However, 
nonperturbatively the situation on the lattice is less clear. Potentially dangerous ultraviolet contributions may appear 
for a naive evaluation of the phase factor of the fermion determinant on the lattice. Despite of these and other potential 
pitfalls, it is our hope that the results of this paper can be compared to lattice QCD and contribute constructively to 
our understanding of QCD at finite density. 
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APPENDIX A: PARAMETERIZATION INDEPENDENCE OF MEAN FIELD RESULTS. 




In this appendix wc illustrate that the mean field result including the 1-loop fluctuations about the saddle point 
(see eq. ^) does not depend on the parametrization of the integration manifold. The nontrivial example we study 
is the phase quenched partition function in the condensed phase, i.e. for ^ > TOtt/S. Because of the presence of an 
exactly massless mode in this phase, we include a pion source which explicitly breaks the U{1) symmetry that is 
responsible for the massless mode. This allows us to use the saddle point approximation to show that the partition 
function is independent of the representation whereas the Jacobian and masses are representation dependent. 

The pion condensate source term enters in the mass matrix 

(Al) 

so that the otherwise massless mode obtains a mass ~ VA. We will consider two different representations of the 
Goldstone fields 

cos Be"' e"f'sme ^ ^ ^^er^^^^,r,^^u. 

\^ -e-"^sin6l cos 6*6-'" J 

where is a vacuum angle that will be determined by the saddle point equations. The Pauli matrices are labeled t^. 
The partition function in either variables is given by (A = XYiV) 

Zi\i.(m,A;/i) = J d6'dad(/)dtJsin6icos6'e2'^''''°'^+2'^™''«™''"™''"'+2^^''^^™'''^™'''^, (A3) 



and 



^ g2m[cos 9 cos tt— sin O^^^^izi] cos aJ+A[cos O^^^^tti +cos tt sin 9] cos uj ( A4) 



where tt = -^/tt^ + 7r| + 7r| . We will evaluate both partition functions in the thermodynamic limit, /t^ — > oo, to — + cx3 
and A ^ oo, where the integrals can be performed by a saddle point approximation. The parameter 9 will be chosen 
such that the saddle point of the tt^ variables is at zero. It turns out that this is also the saddle point of the ^-integral 
in (|A3p . In terms of the 9 — a — (f) — uj variables the squared masses are given by 

rh cos 9, Asin^, 2fL^ sin^ 9 -\ TOCOS0 + Asin^ (A5) 

sin 9 

with 6 determined by the saddle point equation 

cos^sin^ - msin^+ Acos^ = 0. (A6) 

The squared masses in terms TTfc — w- variables can be rewritten by using the saddle point equation which is also given 
by (IMI- We find 

2/i^ + ^^ = ^^, 2/i^sin^ ^+ — TOCOS^ + Asin^. (A7) 

sin 9 cos 9 sin 9 sin 9 

The product of the squared masses is different for the two parameterizations. The contribution from the Jacobian in 
(lASp . sin ^ cos ^ exactly compensates for this difference such that the mean field partition function ([7]) is identical in 
the two cases. In the limit of A ^ to and X fi^ we obtain 

Zi+i.(TO,A;A) = _ ^p:^+^/2p^^ (A8) 

TOY Xsjl^p- — rh? /2(P 

in agreement with the general expression ([7]). 

What we have learned from this example is that the chemical potential dependence of the partition function 
originates both from both the masses and the Jacobian with neither of them being representation independent. 



20 



APPENDIX B: INTEGRALS APPEARING IN CAUCHY TRANSFORMS 



In this appendix we evaluate some integrals that appear in the Cauchy transforms X^''^ given in (I^S]) . 
To simplify the integrals appearing in X^*^^ the following integral is useful 

/■°° 2 2 1 a^+t^ ^ ab\ 

I dtte-P * [Ko{at)Io{bt) + h{at)K^{bt)] - —e^Ko i-^j . 

This identity can be proved by considering the integral 

) POO 

dtte^P^*^[Ko{at)Io{bt) + Io{at)Ko{bt)] = lim / dtte-P^*\K^{at)I^{bt) + I-^{at)K_^{bt)], 



and using the identity 



I sm(^'7r) 



(Bl) 



(B2) 



(B3) 



to replace the K±i^ function by a I±v functions. 

By differentiation with respect to the parameters of this integral, we can derive the following useful identities: 



n 2,2 1 + l 

dtt^e-P * \Ka{at)h{bt) - h(at)KAbt)\ ^ -— 



dtt^e-P^*\Kaiat)lQ{bt) + lQ{at)Kaibt)] ^ 

i 

dtt^e-p''''[Ko{at)Ii{bt) - h[at)Kx{bt)] = 



4p4 



bK. 



(2 



ab 
V 



aKi 



)Ko 



ab 
V 



ab 



a(a'^ + 35^) 
V 



afe / ab 



a6 



(B4) 
(B5) 

(B6) 



APPENDIX C: ASYMPTOTIC EXPANSION OF INTEGRALS OCCURRING IN THE X'*' 

In this section we calculate the asymptotic expansion of the integrals occurring in the X'*^) given in eq. (|46|) . We 
separately consider the integral 



dxa^Pexp[-^]ifo [j^ ] 



(CI) 



in the normal phase and in the condensed phase. 

In the normal phase, for m/4 < ji^ < rh/2, the saddle point at a; = 4/i^ — m of the integral (jCip is inside the 
integration domain. For m ^ oo and fi^ oo the Bessel functions can be expanded to leading order resulting in 



^normal 



2m 



(C2) 



For fi^ < m/4 the leading saddle point contribution comes from the second integral in the X'-''\ We leave it up to 
the reader to show that the leading order saddle point approximation results in the same expression. 

In the condensed phase the leading order expansion of Sp^i, cancels in the expression for the average phase. Therefore 
we have to include the subleading terms in its asymptotic expansion. Including the expansion of the Bessel functions 
to subleading order we obtain for Sp,i, 



s: 



condensed 



■'/Sfi -xm/ifi-'+x 



2X771 



(C3) 



Next put X ~ 771 — t. Since the integral is dominated by the vicinity of f = we can safely extend the integration 
range to oo, 
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To leading nonvanishing order it suffices to expand exp(— as 1 
find 



8^ 



After performing the integral over t we 



s: 



condensed 



-e m' 



2m 

2 - 4S 



2m 



4(P-1) 

(2-4)2 ji^{2~m/ji^f (2-m//i2)^m 



(C5) 



Combining the different terms we obtain for the thermodynamic limit of the expectation value of the quenched phase 



(exp(2i6'))jV/=o 



1 1 / 2m-m/2A^-2/i^ 



\/27rm A^(2 - m/fi^Y V m 



for 



1 < rn < 2/i 



(C6) 



The condition for the validity of the derivation of (jC6P is that (2/1^ — m)//i ^ 1, which is only violated very close to 
the critical point because the natural magnitude of this ratio of 0{^/V). We can also evaluate the integral (jC4p when 
this condition is not satisfied. Then the integral (jC4p is cut-off by the the Gaussian factor exp(— i^/2/x^) instead of 
the exponential factor resulting in an expression that remains finite for m = 2/i^. 
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